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Abstract 

We investigate the effect on the sphaleron in the two Higgs doubiet 
eiectroweak theory of including CP violation in the Higgs potential. To 
have better control over the relation between the sphaleron energy and the 
physical quantities in the theory, we show how to parametrize the Higgs 
potential in terms of physical masses and mixing angles, one of which causes 
CP violation. By altering this CP violating angle (and keeping the other 
physical quantities fixed) the sphaleron energy increases by up to 10%. We 
also calculate the static minimum energy path between adjacent vacua as a 
function of Chern-Simons number, using the method of gradient flow. The 
only effect CP violation has on the barrier is the change in height. As a by- 
product of our work on parametrization of the potential, we demonstrate 
that CP violation in the Higgs sector favours nearly degenerate light Higgs 
masses. 



1 Introduction 

There has been much recent interest in the possibility of generating the asym- 
metry in baryon number observed today during an eiectroweak phase transition 
[|l], [2j. For a theory to allow any baryon asymmetry generation, it must satisfy 
certain conditions, originally identified by Sakharov ||: a sufficient departure 
from thermal equilibrium, violation of charge conjugation (C) invariance and the 
combination of C with parity (P) invariance, and violation of baryon number 
conservation (B). 
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The standard model of Weinberg and Salam posseses all three properties. C 
is violated maximally in that only left handed fermions transform nontrivially 
under SU(2). CP is violated by a small amount given by the phase in the CKM 
matrix, and so CP is violated only through the Yukawa couplings. B is conserved 
to all orders in perturbation theory at zero temperature, but can be violated 
non-perturbatively either through the quantum tunneling of instantons between 
inequivalent vacua [|] or through the formation and decay of sphalerons || at 
finite temperature ||. The required departure from thermal equilibrium occurs 
during the sufficiently first order phase transition that breaks the gauge SUl(2) x 

Uy(l) ->U em (l). 

Unfortunately, although the SM has all the ingredients for baryon asymmetry 
generation, it does not seem able to create the observed B asymmetry; the CP 
violation is too small, and the phase transition is not sufficiently first order, and 
may even not be present at all (see for example J7| and references therein). 

Attention has turned to extentions of the standard model. Of particular 
interest has been the two Higgs standard model, 2HSM. This is a well motivated 
extension in that it allows sources of CP violation to enter through the Higgs 
potential, and it includes the Minimal Supersymmetric Standard Model as a 
subset. 

It is the 2HSM model that we consider throughout. Curiously, although 
the sphaleron has been studied before in the 2HSM ||, [9|, [TO], the effect of CP 
violation in the Higgs potential has not been studied in detail. The purpose of 
this paper is to investigate how the sphaleron reacts to CP violation in the Higgs 
sector. In doing so we must decide how to parametrize CP violation, as this is 
by no means unique. The most satisfactory solution, in our opinion, is the most 
physical one: we parametrize CP violation by the mixing angle between the CP- 
odd and CP-even Higgs states (by CP-even and CP-odd we mean the quantum 
numbers the states would possess in the absence of mixing). This leads us to 
the question of how we fix all 10 independent parameters in the most general 
SU(2) xU(l)-invariant renormalisable two Higgs potential with acceptable levels 
of flavour- changing neutral currents. We solve this problem by inferring them 
from the Higgs masses and mixing angles, and the Higgs vacuum expectation 
value. This leaves 2 parameters we must fix by hand, one CP violating, and the 
other a Higgs 4 point coupling. This way of proceeding makes it much easier to 
satisfy experimental constraints and vacuum stability from the outset, although 
we must check that the resulting parameters also satisfy the constraints from 
boundedness of the potential. 

Intriguingly, we find that it is quite difficult to cover a wide range of the CP- 
violating mixing angle, unless the CP-even states are almost degenerate in mass 
and the CP-odd state is very massive . However, once we choose such a set of 
masses (100, 110, and 500 GeV), the sphaleron energy varies by as much as 10% 
as the CP-odd mixing angle changes from to about tt/3. This has important 
implications for the theory of electroweak baryogenesis. There are strong up- 
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per bounds on the lightest Higgs mass, which derive from the requirement that 
sphaleron transitions not wipe out any generated baryon asymmetry in the bro- 
ken phase, which in turn constrains the strength of the transition. Even in the 
MSSM these bounds are close to being violated by the experimental lower bound 
on the Higgs mass |TT . 



Another novelty of our approach is the use of the gradient flow method [12] to 



find the entire barrier, that is we calculate numerically the static minimum energy 
path (SMEP) between adjacent vacua as a function of the Chern-Simons number. 
This has not been done before for the two Higgs model, and we were motivated 
to do this to search for signs of the underlying CP asymmetry in the barrier. 
The gradient flow technique is slightly easier to implement than methods which 
calculate the energy subject to a constraint |13|, [HJ, and gives an intuitive feel 
for how the field might relax to the vacuum after a sphaleron has been created. 

In Section 2, we introduce the model and our parameter definitions for the 
two Higgs potential. We restrict our attention to SU(2), as we then need consider 
only spherically symmetric field configurations. This is sufficient to demonstrate 
the point that one should consider more realistic Higgs potentials than hitherto. 
In Section 3 we construct a general spherically symmetric ansatz for the sphaleron 
which allows for CP-violation. In Section 4 we show how to determine the cou- 
plings of our potential as a function of the physical masses and mixing angles, 
and derive the bounds on these couplings from vacuum stability and boundedness 
of potential. We verify the sphaleron solution for a general 2HSM in Section 5 
and proceed to find the entire barrier between neighbouring vacua by gradient 
flow. We then find the height of the energy barrier as a function of the only 
CP-violating physical parameter, 9, the mixing between CP-even and CP-odd 
neutral Higgses. In Section 7 we supply our conclusions and outline possibilities 
for future work. 



2 The Model 

We consider the bosonic sector of the electroweak lagrangian with two Higgs 
doublets in the limit 9 W = 0, so that 

c = -\f; u f^ + (D,<p a y(Dy a ) - vfa, <p 2 ) (1) 

where 

F% = W-^; + ^X c (2) 

D^ a = (dn + gWJta (3) 

W, = Wf-^ (4) 

Higgs fields are labelled by a = 1,2. The most general renormalizable potential 
for two Higgs doublets has 14 real couplings. To reduce the effect of flavour 
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changing neutral currents the symmetry <p\ 



b 2 is often imposed 



but relaxed for dimension two terms fll5| . Such a potential will have ten real 
couplings, and can be written 



Vtyufa) = (A1 + A3) +(A 2 + A 3 ) [4 
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Vg 

V2 



J8 a 



(5) 



(6) 



where Q\ = and 9 2 = £. We are free to redefine our fields, and we chose to 
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The potential can now be written as a function of nine real 
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The tenth coupling has been shifted in the redefinition to the CP transforma- 
tion properties of the fields, so that now under a CP transformation <p a — > 

We chose Eq. |7| as our potential and set £ = throughout, so that our vacuum 
configuration is invariant under CP, and the only source of CP violation is in the 
coupling X2- The physical effect of \2 will be to produce a non-zero mixing angle 
between the CP-even and CP-odd neutral Higgses. 



3 The Sphaleron Ansatz 

A sphaleron is a static unstable finite energy field configuration, which represents 
the top of a static minimum energy path, SMEP, connecting inequivalent vacua. 
Each vacuum is distinguished in the unitary gauge by its Chern-Simons number 
ncs, defined as 



- ^ f A 1„TfO 



1 



WfdjWg + -ge ahc W*W]Wl 



(12) 



32tt j 



d 6 xK\ (13) 



where d^K^ 1 = F^ u F atlu . The vacua have integer rics; and are related by a large 
gauge transformation. 

In the standard model there are twelve U(l) currents one for each left- 
handed fermionic species, defined as J® = \l/^7 M \I / ^' ) . At the quantum level, 
these currents are not conserved: 

32vr 2 ^ v 1 

Hence the fermion numbers N^ 1 ' are not conserved, by an amount which by the 
change in the Chern-Simons number. We can infer then that the baryon and 
lepton number violation as the background gauge fields change is 

AB = AL = 3An cs (15) 

Vacua with different Chern-Simons number are connected in configuration space 
along a path whose energy is always finite, so it is possible, with sufficient energy, 
to change B+L classically. The least energetic way of doing this is to form a 
sphaleron. 

Motivated by |J, we chose a more general version of the spherically symmetric 
ansatz of [ TH j : 

<P a = ^(F a + tG a x a a a ) (16) 

W = -A x a ^ (17) 
g 2% 
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Wi = - 
9 



[1 + 13) . a . M „ 

^ £-aim&m + ^ Kpai %a%i) + A\X a X\ 



2i 



18) 



where F Q = a Q + and G a = c a + id a . 

On substituting our ansatz into Eq. [I], our theory posseses spherical symmetry 
only if 

F a = A(r, t)G a (19) 

we assumed spherical symmetry, and verified that all terms that were zero because 
of Eq. [19] did indeed vanish. With this condition the A4 term disappears from the 
potential. The A4 term is solely responsible for the mass of the physical charged 
Higgses, rriH± = X4V 2 /2. For our ansatz, imposing spherical symmetry forces the 
A4 term in the potential to zero. Thus the charged Higgses decouple, and the 
choice of A4 is irrelevant for the sphaleron. 



4 Choosing Couplings 

We expanded 

^ = 7=2 



v a + H c n 



(20) 



to get the mass-squared matrix M G in the gauge basis he of neutral scalar Higgses, 
where 

h G = [ReH°, Re# 2 °, A ] . (21) 

Here, A = (ImH® cos /3— Imi?° sin (3), and is odd under CP, the other two neutral 
Higgses being even. The components of the mass matrix in the gauge basis are 



M G (1,1) 




4(Ai + X 3 )cos 2 (3+ (A 4 


- + Xi) sin 2 /? 


V 

Y 


(22) 


M G (1,2) 


= M G (2, 1) = (4A 3 + A + 4 


- xi) cos P sm 


4 


(23) 


M G (1,3) 


v 2 

= M G (3,1) = X 2smP- 






(24) 


M G (2,2) 




4(A 2 + A 3 )sin 2 /3 + (A + 


+ Xi)cos 2 /3 


v 2 
2 


(25) 


M G (2,3) 


v 2 

= M G (3,2) = X2COS p- 






(26) 


M G (3,3) 


v 2 

= (A + -Xi)y 






(27) 



We denote the neutral Higgs fields and the mass-squared matrix in the physical 
basis by hp and M P respectively, with hp = [h 2 , hi, /i 3 ]and M P = diag [m 2 , m 2 , m|]. 
The reason for our ordering convention is that in the absence of CP violation the 
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lightest Higgs is conventionally the second entry in the physical higgs state vector 
hp. We chose as our three mixing angles the Euler angles if), 9, and 0, such that 



hp 
M P 



R z (il>)R y (6)R z (4>) h G = Dty, 9, 
D(^, 9, 0) M G ZT 1 ^, 9, 0) 



'G 



(28) 
(29) 



where 



R*(a) 



cos a sin a 
— sin a cos a 
1 



R y (a) 



cos a — sin a 

1 
sin a cos a 



(30) 



We see that; is responsible for the mixing between the CP even Higgses to give 
H and h°, (H° = ReH° cos + ReH$ sin 0, and h° = -ReflJ sin + Re# 2 ° cos 0) . 
The angle is responsible entirely for the CP violation as it mixes the CP-odd ^4° 
with the CP-even H°, and ip then shares this CP violation between H° and h°. 
For the case 8 = 0, there is no CP violation; mi = rriho, m<i = m^o, iris = m^o, 
and the mixing angle between the two CP-even states is + ip. 

If we invert Eq. we can in principle find Mq as a function of the physical 
parameters, which we call v 2 X(mi,m2,m3,<f),8,if}). The matrix Mq is itself a 
function of eight unknown parameters (Ai, A 2 , A 3 , A 4 , A + , Xi, X2, P) and one known 
parameter (v = 246 GeV), where 



v 



v\ + vl, 



tan (3 



V2 



(31) 



The procedure is to find as many as posible of the parameters (Ai, A2, A3, A4, A + , xi, 
X21 /3) in terms of the physical quantities, which are four masses and three mixing 
angles (m 1; m 2 , m 3 , mji±, 0, 9, ip). Clearly, we have to fix one of them by hand: 
this we chose to be A3. 

The charged Higgs mass m#± is determined only by v and A4, and will only 
be relevant when we consider the boundedness of the potential and the stability 
of the vacuum. We thus now have seven unknown couplings Ai, A 2 , A3, Xu X2, P- 
From Eq. 123 and Eq. EM 



X2 

(3 = arctan 



2 V /X(1,3) 2 + X(2,3) 2 , 



X(2,3) 



(32) 
(33) 



In the limit 9 = 0; X(l, 3) = X(2,3) = 0, X2 = 0, and (3 can be taken as the 
limiting value of the ratio. All other couplings are determined as below. 

Using the non-independent relation Tr Mq = m\ + m?, + m 2 , and eliminating 
A + using Eq. ^7], we were able to solve Eqs. [22|, ^3| and |5| to determine all the 
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couplings. On the understanding that A3 has already been fixed by hand, we get 

1 



Ax 
A 2 
A + 

Xi 



[(X(l, 1) cos (3 - X(l, 2) sin (3) - 2A 3 cos 2(3 sin (3} 
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-2A, - 



X(3,3)+X(l,2)- 
X(3,3) + X(l,2 



sin /3 cos (3 
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(34) 
(35) 
(36) 
(37) 



sin (3 cos /3 

These are linear equations which are easily solved. 

We now turn to the constaints that the parameters so derived must satisfy. 
There are eight conditions on our potential, which derive from its boundedness 
and the stability of the vacuum state. For boundedness of the potential we require 
that the eigenvalues of 



dxidxj 

are all positive, where i,j = 1,2,3,4; x± = 4> 
£4 = Im(0j0 2 )- This gives the conditions 



101, X 2 



202, X 3 
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For a stable vacuum, we require 



> xi + xi 



ml > 0, m\> 0, 



mi > 0, 



m H ± > 0. 



(38) 
and 

(39) 
(40) 
(41) 
(42) 

(43) 



The last of these conditions is easily satisfied by A 4 > 0. 

On substituting Eqs. |32|-^7| into the inequalities [39|-fH we could derive con- 
ditions directly on masses and mixing angles. In practice, we picked masses and 
mixings, calculated couplings for a suitable A 3 , and then verified that [39|-442| held. 
Although we always used values of masses and mixing angles that satisfied all the 
conditions for boundedness and stability, in general it was quite difficult to find 



masses and mixings such that the inequalities |39|-442| were satisfied. For example 
for mi=100 GeV, m 2 =300 GeV, m 3 =400 GeV, m H ±=50 GeV, V=0=f , the only 
allowed range of 9 is 2.40 < 9 < 2.85. 



5 The Barrier 

We wanted to find the SMEP between vacua differing in n cs by 1. This has been 
done in the one doublet case using the method of Lagrange multipliers to fix a 
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constraint, either ncs itself ||13||, or the distance in field space from the sphaleron 



14] . We chose instead gradient flow [ 12| . Gradient flow is defined from the static 



energy £ s by 

m = ~ K JJ (44) 

where k is a friction term, and / are the fields of the theory. In our case / = 
(a a , b a , c a , d a , a, ft, A\) (note that in the static energy, A = 0). Gradient flow 
will always describe a path of minimum energy since it forces the fields to evolve 
with a velocity orthogonal to the contours of £ s . 

To flow down the barrier we first need to find the sphaleron. We chose to 
work in the radial gauge, A\ = 0, which is the most convenient, for in this gauge 
so fixing the boundary conditions of all fields fixes ncs- To start at the sphaleron, 
we fixed the boundary conditions corresponding to ncs — \i which are 

fla(0) = 6 Q (0) = c Q (0) = rf Q = a(0) = 0, 0(0) = 1, (45) 

a a (oo) = b a (oo) = d a (oo) = a(oo) = 0,(00) = ft (00) = 1 (46) 

We used Simultaneous Over Relaxation, and Chebyshev Acceleration (without 
even-odd ordering) |17| to relax to the minimum energy configuration, and found 
a solution consistent with || for the same values of the parameters, in that we 
had the same energy and field profiles were indistinguishable when compared by 
eye. We used a grid size of 201, and considered our fields to have converged to 
solutions when the change in the absolute value of the fields integrated over the 
grid was < 10~ 13 . 

This gave us the initial condition for our gradient flow. A technical problem 
now arises, because once the fields start flowing, A x has an equation of motion 
and will not in general remain zero, that is, the configuration will not remain in 
the radial gauge. It is useful to stay in the radial gauge in order to be able to 
compute the Chern-Simons number easily. Hence, in order to keep in the radial 
gauge we carried out a gauge transformation 

a a + ic a = e~ l t (a a + ic a ) (47) 

b a + id a = e~ l %(b a + id a ) (48) 

a + ift a = e- ie (a + ift) (49) 

A x = Ai-Q' (50) 

after each step of evolution of gradient flow, where we chose so that 

= / A x dx (51) 



Naturally, there was no evolution unless we put a small perturbation in one of the 
fields. For a small positive perturbation we arrived at a vacuum with ncs = 1, 
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for a small negative perturbation we arrived at a vacuum with ncs — 0. The 
barrier was independent of choices of k, and we used k = 0.2 throughout. 

We checked our code by computing the barrier and the field profiles for the 
sphaleron in the one Higgs model, and comparing by eye the results of Nolte and 
Kunz jrjL finding no noticeable differences. 
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Figure 1: The sphaleron barrier for mi=100 GeV, m 2 =300 GeV, m 3 =500 GeV, 
if)=j, > 0=-f and A3=0. The solid, dot-dashed and dashed lines are the 
contributions to the energy from the gauge field, the mixed or covariant derivative, 
and potential terms respectively. 



The result, for one particular choice of masses and mixing angles, is shown 
in Fig. [I]. The barrier showed no unusual features for having an extra doublet, 
and CP violation. In particular, it is symmetric around Chern-Simons number 
|. That this should be so is not obvious a priori, because ncs is odd under CP, 
and moreover changes by integer multiples under gauge transformations. Hence, 
a combination of a CP transformation and a large gauge transformation changes 
ncs to 1 — rics- It is therefore interesting that the barrier should be symmetric 
under this interchange when the underlying theory is not. It is also interesting 
to note the small contribution of the potential energy. 



6 Energy as a function of CP violating angle 

We were unable to find sensible values of mi, m 2 , m%, mjj±, 4> and ip, that satisfied 
Eqs. |39|-f42| for the entire range < 9 < tt. In turned out that the largest ranges 



of 9 favoured degenerate mi and 777,2, large 7713, and small m H ±. We plotted height 
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of the energy barrier as a function of 9 for the allowed range. In Fig. |2| we display 
the energy as a function of the CP violating angle 9, for m^lOO GeV, m 2 =110 
GeV, m3=500 GeV, ip — |, = | and A3=0, values which are consistent with 
current bounds on Higgs masses |TT 
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Figure 2: The sphaleron energy as a function of the CP violating Higgs mixing 
angle 9, for Higgs masses m^lOO GeV, m 2 =110 GeV, m 3 =500 GeV, and CP- 
conserving mixing angles tp = ~, <p = |. The only coupling undetermined by the 
masses and the mixing angles is A3, which was set to zero. The solid, dot-dashed 
and dashed lines are the contributions to the energy from the gauge field, the 
mixed or covariant derivative, and potential terms respectively. 

We note that the energy changes from 9 to 10 TeV over the range of 9 we could 
obtain, or by roughly 10%. This is not insignificant, in view of the tightness of the 
upper bound on the lightest Higgs mass, which comes from the requirement that 
the sphaleron rate be sufficiently low in the symmetry-broken phase at the bubble 
nucleation temperature T&. Increasing the sphaleron energy by this amount would 
mean that the ratio v(Tb)/Tb can be smaller, and hence that the transition can 
be more weakly first order. This means that the lightest Higgs is allowed to be 
more massive (see, e.g., U for a more detailed explanation). 

Clearly, a more precise determination of the bounds must take into account 
finite temperature corrections, and use the true SU(2)xU(l) theory. At the same 
time one should compute the ratio v{T b )/T b with the correct Higgs parameters. 
It is interesting to note that existing lattice calculations assume that the 
Higgses in two-doublet models are not degenerate in mass, in order that all but 
the lightest can be integrated out. We find that in order to get large range 
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of CP-violating angles, one must have almost degenerate light Higges. It may 
well be worth recalculating the strength of the transition as a function of the 
(zero-temperature) Higgs masses in the almost degenerate case. 

7 Conclusions 

In this paper we have made a start on investigating the effect of CP violation in 
the Higgs potential on the sphaleron energy, by looking at the sphaleron in the 
pure SU(2) two Higgs theory. We found firstly that the sphaleron energy changed 
by about 10% as we scanned though a range of our CP- violating parameter, which 
we defined to be the mixing of the CP-odd with the CP-even Higgs states. This 
indicates that it is worth going on to calculate the sphaleron energy with finite 
temperature corrections taken into account, and in the full SU(2) xU(l) two Higgs 
theory. In this way the bounds deriving from the requirement that the baryon 
asymmetry not be destroyed in the symmetry-broken phase can be more precisely 
worked out. We also used a new method, that of gradient flow, to follow the fields 
as the sphaleron decays to the vacuum, and to exhibit the energy as a function of 
Chern-Simons number. The shape of the barrier shows no qualitative difference 
from that in CP-conserving theories: in particular, it is symmetric around Chern- 
Simons number |. In a CP- violating theory, this feature is not a priori obvious, 
and perhaps points to a deeper symmetry in the equations. 

Finally, in order to relate our sphaleron energies to physical quantities, we 
adopted a novel strategy of computing as many as possible of the parameters 
of the Higgs potential from the masses and mixing angles of the Higgses. With 
four masses, three mixing angles, and the Higgs vacuum expectation value we 
needed to choose only two of the parameters without physical input. We found 
as a by-product that a large range of the CP-violating parameter could only be 
obtained if two of the Higgs were nearly degenerate in mass, and one much more 
massive. This again has important implications for the lattice calculations of the 
strength of the electroweak phase transition, which hitherto have assumed that 
the Higgses are well separated in mass. 
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